Phases intermediate between the two dimensional electron liquid and Wigner crystal 
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We show that there can be no direct first order transition between a Fermi liquid and an in- 
sulating electronic (Wigner) crystalline phase in a clean two-dimensional electron gas in a metal- 
oxide-semiconductor field-effect transistor (MOSFET); rather, there must always exist intermediate 
"micro-emulsion" phases, and an accompanying sequence of continuous phase transitions. Among 
the intermediate phases which we find are a variety of electronic liquid crystalline phases, including 
stripe-related analogues of classical smectics and nematics. The existence of these phases can be 
established in the neighborhood of the phase boundaries on the basis of an asymptotically exact anal- 
ysis, and reasonable estimates can be made concerning the ranges of electron densities and device 
geometries in which they exist. They likely occur in clean Si MOSFETs in the range of densities in 
which an "apparent metal to insulator transition" has been observed in existing experiments. We 
also point out that, in analogy with the Pomaranchuk effect in 3 He, the Wigner crystalline phase 
has higher spin entropy than the liquid phase, leading to an increasing tendency to crystallization 
with increasing temperature! 



In discussions of the theory of the two dimensional elec- 
tron gas (2DEG), it is generally accepted that, as a func- 
tion of electron density n, there is a first order quantum 
(T = 0) phase transition from a high density liquidQ] 
to a low density Wigner crystalline phase0- This as- 
sumption is reasonable in the case of a triangular Wigner 
crystal due to the presence of cubic invariants in the Lan- 
dau free energy jS, and for other lattices due to the gen- 
eral expectation 6] that fluctuations will always render a 
freezing transition first order. The transition is thought 
to occur when the dimcnsionless ratio r s = [^(ob) 2 ] -1 / 2 
exceeds a critical value r s = r c ~ 38, where as is the 
effective Bohr radius in the semiconductor. However, this 
generally accepted picture is manifestly incorrect for the 
2DEG in a metal-oxide-semiconductor field effect tran- 
sistor (MOSFET), and possibly more generally! 

Each electron in the 2DEG in a clean MOSFET drags 
along with it an image charge in the ground-plane above. 
Consequently, at small separations, the interaction be- 
tween the electrons is the V(r) ~ e 2 /er Coulomb inter- 
action, while for separations larger than the distance to 
the gate, d, it is the repulsive dipole-dipole interaction, 
V(r) ~ 4e 2 d 2 /er 3 . (Here e is the dielectric constant of 
the host semiconductor.) In 2D systems with dipolar 
interactions, the following simple argument leads to the 
concussion that first order phase transitions are forbid- 
den: In systems with interactions that fall more rapidly 
than 1/r 2 , there exists a "forbidden" range of densities in 
the neighborhood of a first order phase transition where 
macroscopic phase separation reduces the free energy of 
the system. However, when we come to compute the 
surface tension between two macroscopic phases, we find 
that 1/r 3 interactions are marginal: for shorter range in- 
teractions, there is a well defined scale independent sur- 



face tension, a, while for longer range interactions, a is 
scale dependent. Specifically, for dipolar interactions, the 
interfacial contribution to the free-energy of an arbitrary 
macroscopic mixture of two phases is (see, e.g. 0,0) 
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Here, the arclength integral, ds, runs along 8] the inter- 
faces between the two phases, 0(s) is the local orientation 
of the interface, (Tq(6) is the (in general orientation de- 
pendent and by assumption positive) short-range piece 
of the surface tension, dl runs along the interfaces @, 
crj = 2e 2 (An) 2 G? 2 /e, and An is the density difference 
between the coexisting phases. The second (non-local) 
term in Eq. ^ comes from the long-range parts of the 
dipolar interaction. One can also view it as the leading 
finite size correction to the capacitance of parallel-plate 
capacitors due to the fringing fields [9j. 

It is important to note that the second term in Eq. ^ 
gives a negative contribution to the effective surface ten- 
sion which diverges logarithmically with length; for ex- 
ample, an isolated straight segment of interface of length 
L has F<j = L {cto — <j\ log[L/2c?]}. This implies that 
there is an absolute instability of the macroscopically 
phase separated state - in the regime of the phase dia- 
gram where a classical Maxwell construction would lead 
to two-phase coexistence, a state formed from a "mi- 
croemulsion" of the two phases (with a character and 
length scale to be determined), has lower free-energy! 
Thus, instead of a first order transition between two 
phases, there must always be an intermediate regime in 
which one or more microemulsion phase occurs, bounded 
by one or more line of continuous phase transitions. 

At this point we would like to compare this situation 
with the Coulomb case (no ground plane) where macro- 
scopic phase separation is forbidden. The nature of the 
phases that result from the "Coulomb frustrated phase 
separation^! in what would otherwise have been 
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the forbidden range of densities is an issue of potentially 
relevance in many highly correlated materials. However, 
the inhomogeneities that occur in this situation are typ- 
ically microscopic in scale, and so difficult to distinguish 
from more familiar charge density wave structures ^ij- 
Moreover, the relevant microscopic details are difficult to 
treat with any degree of rigor, (ft is an interesting^^ 
question, which we would like to reopen, whether there 
are intermediate phases between the Fermi liquid and 
Wigner crystal phases in the 2DEG with pure Coulomb 
interactions.) 

The character of the microemuslion of the two coexist- 
ing phases is determined by minimizing F a in Eq. ^ the 
result depends on how anisotropic the function ao(9) is. 
The case where cto (9) is independent of 9 has been consid- 
ered in different contexts, including lipid films (e.g. Ref. 
|l2j1. two dimensional uniaxial ferromagnets (e.g. |13|1. 
and the 2DEG in MOSFET's [3 [T^. The resulting 
phase diagram includes both stripe and bubble phases, 
with stripes preferred in the center of the phase separated 
region and bubbles generally thought to be slightly lower 
in energy when one phase is in extreme minority. Cur- 
rent estimates [5| place the difference between the dilute 
stripe and bubble energies at about 6%. In the earlier 
literature, it was assumed [l2l IT^| that there is a direct 
first order transition between uniform stripe and bubble 
phases. This is incorrect, even at mean-field level, since, 
as we have shown, first order phase transitions are for- 
bidden. Thus, a sequence of continuous phase transitions 
(which we discuss below) must replace the putative first 
order transition |14| . 

In the present case, where at least one of the two co- 
existing phases is crystalline, the angular dependence of 
<To(9) is not negligible, reflecting the tendency of crys- 
tals to facet. Clearly, a strong angle dependence of <Jq(9) 
tends to favor stripe phases (where all interfaces lie along 
the direction in which gq(6) is minimal) relative to any 
form of bubble phase. 

In the present paper, we characterize the phase dia- 
gram, and in particular the universal aspects of the inter- 
mediate phases and phase transitions that are expected 
at low or zero temperature in an ideal MOSFET (i.e. 
in the absence any disorder). We will consider explicitly 
the case in which d is large compared to the spacing be- 
tween electrons, nd 2 ^S> 1, as in this limit (as we shall 
see) fluctuation effects are parametrically small and an 
appropriate mean-field theory provides a valid zeroth- or- 
der description of the phases. In Sec. [fl we first discuss 
the mean-field phase diagram, then in Sees. Illland lTTll we 
discuss the effects of weak thermal and quantum fluctu- 
ations, respectively. In Sec. IIVI we discuss some of the 
implications of the present results for the properties of 
real devices (which, alas, have non-negligible disorder), 
and in Sec. [V] we discuss some incompletely developed 
ideas concerning further implications of the present line 
of analysis. 



I. MEAN-FIELD PHASE DIAGRAM 

Two dimensionless parameters determine the physics 
of the 2DEG in a MOSFET, r s (defined above) and a B /d. 
Let us start with a discussion of the zero temperature 
mean-field phase diagram of this system, assuming only 
uniform states. If nd 2 ^> 1, the free-energy per unit area 
can be represented by the sum f(n) = f {C) +f {el) of the 
energy density of a capacitor = (en) 2 /2C and the 
internal free-energy density of the electron liquid f( el \ 
Here C = (eeQ -1 is the capacitance per unit area. At 
high electron densities, r s « 1, the kinetic energy of 
the electrons is much larger than their potential energy, 
so the system forms a Fermi liquid. At small densities 
r s ^S> 1 (but still nd 2 3> I ) the Coulomb energy of the 
electrons is much larger than the kinetic energy, so the 
ground state is crystalline. 

However, at even smaller densities when nd 2 <C I, the 
electrons interact only via dipole interactions, so the ki- 
netic energy is larger than the potential, and the sys- 
tem again has a Fermi liquid groundstate. (See dis- 
cussion surrounding Eq. 1121 1 For d/as S> 1, this im- 
plies that the phase diagram of the system has reen- 
trant transitions as a function of n (along the dashed- 
dotted trajectory in Fig. 1) from a Fermi liquid phase 
for n > n c r~ 2 (7ra^) _1 to a Wigner crystal phase 
for n c > n > n c \ ~ (ndr) , to a Fermi liquid phase 
for n c \ > n. With decreasing d/as, n c i and n c move 
toward each other, until for d < d c ~ r c aB, the Wigner- 
crystal phase disappears entirely. This is represented by 
the dashed line in Fig. 1. 

As a next step, we improve this phase diagram by al- 
lowing for the possibility of inhomogeneous states. There 
is a range of forbidden densities about the critical den- 
sity in which macroscopic phase separation into regions 
of high and low density phase has lower free energy than 
the uniform state. 

Let us briefly review the salient features of the Maxwell 
construction for phase coexistence, as applied in the 
present context. For given average density, n, we con- 
sider a state in which a fraction, x, of the system is at 
a higher than average density, n+ > n, and a fraction 
(1 — x) is at a lower than average density, n_ < n, such 
that xn + + (1 — x)n- = n. We then minimize the total 
free energy with respect to n + and n_ . The result of this 
minimization is an implicit expression for the densities of 
the two coexisting phases, 

n+ _ n- _ [/(n_) - f(n+)} 
^ + + -C-^- + —- An ' (2) 

where fi± — —df( el >(n±)/dn± are the chemical poten- 
tials in the two phases, and An = [n + — nJ\. Phase 
coexistence occurs for n_ < n < n +1 where the fraction 
of the two phases is determined by the lever rule, 

x = (n — n_)/An. (3) 

Eq. [21 is somewhat complicated, but it can be greatly 
simplified when the forbidden region is relatively small 
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FIG. 1: The T = phase diagram of the 2DEG in an MOS- 
FET. The dashed line indicates the mean-field critical density, 
n c (d), where the free energies of the uniform Wigner crystal 
(W.C.) and Fermi liquid (F.L.) phases cross. The solid lines 
mark the boundaries of the regime of the intermediate mi- 
croemulsion (stripe or bubble) phases. At mean-field level, 
these solid lines are Lifshitz transitions. They approximately 
coincide with the regime of macroscopic two-phase coexis- 
tence (n_ < n < n+) derived from a Maxwell construction. 
The hatched area represents the regime in which the regions 
of the two coexisting phases have sizes of order the electron 
spacing, so quantum fluctuations are order 1, and hence may 
substantially alter the mean-field character of the phases and 
phase transitions. The cross-hatched areas denote the regimes 
of Coulomb frustrated phase separation where even the mean- 
field character of the phase diagram is not known. 

(An <C n c )', in this case, we can linearize the density 
dependence of the free energy about the critical density, 

& l) (n±) = / (eZ) K) - M± (n± - O + . . . (4) 

where . . . represents higher order terms in powers of 
(n± — n c ). To this level of approximation, 

n ± =n c ±^; An = (/ *~ "^K (5) 

The discontinuity of the chemical potential, (//_ — fi+) > 
0, is determined by microscopic physics, and is only small 
to the extent that the putative transition is weakly first 
order. Whether or not the transition is strongly first 
order, for d large, An is self-consistently small. 

The validity of the Maxwell construction rests on the 
implicit assumption that the interface energy between 
the coexisting phases is positive, so the amount of inter- 
face is minimized. As we have seen, in the dipolar case 
this assumption is invalid. We can construct a state with 
lower free energy by making an inhomogeneous mixture 
of the two coexisting phases to increase the amount of in- 
terface. To complete the mean field analysis, one should 
minimize Eq. ^with respect to the shape of the minority 
phase regions at given area of the phase - the area being 
given, to first approximation, by the Maxwell rule. 



A. Stripe Phases 

To begin with, let us consider only striped phases. This 
is fully justified in the case of strong anisotropy of the 
surface energy. (As we will see letter even in the oppo- 
site case when a (6) is isotropic, there are regions in the 
phase diagram where this assumption is relevant.) The 
interfacial free energy density for a striped phase is easily 
computed from Eq. [2 to be 

f a = L~ 1 {2(70-4(71 ln[L+L_/dL}}, (6) 

where L± are the widths of the high and low density 
regions, respectively, and L = L+ + L_ is the period of 
the stripe structure. Minimizing Eq. at fixed areal 
fraction of the high density phase, x = L+/L, we get 

L-^-e^; L + = *W+7 ( 7) 
x (1 — x) 

with 7 = (To/2(7i. It is important to note that as x — > 0, 
the stripes of the high density phase approach a finite 
limiting width, L+ — > Lq = de 1+7 , although the spacing 
between stripes, L_, diverges in proportion to 1/x. Also, 
because the minimized value of f a — —Aai/L is negative, 
the region of stability of the striped phase in fact extends 
somewhat beyond the edges (n_ and n+) of the two- 
phase region derived from the Maxwell construction. 

Finally, it is necessary to estimate the magnitude of 7; 
if it is of order 1, then Lq ~ d, but if 7 ^ 1, then Lq 
is exponentially larger than atomic lengths. So long as 
the stripe phase occurs in a relatively narrow range of n, 
we can use Eq. [S] to estimate erj , with the result that 
7 ~ aoe 2 /e[ii+ — /i-] 2 , which is a ratio of microscopic 
electronic energies. Thus, except under special circum- 
stances, we expect that 7 ~ 1, and hence that Lq ~ d. 
However, so long as nd 2 ^> 1, the stripe widths are still 
large compared to the spacing between electrons, which 
validates the macroscopic approach taken here. 

In short, at mean field level, as a function of decreasing 
density the system evolves from the Fermi liquid phase, 
through intermediate stripe phases, to the Wigner crys- 
tal, as summarized in Fig. 2a: 

• 1) Starting in the uniform Fermi liquid phase, as 
the density is varied across n+, the system under- 
goes a transition to a stripe phase, consisting of a 
periodic array of far separated stripes of Wigner 
crystal, with characteristic width Lq. This transi- 
tion is analogous to a Lifshitz transition, in that the 
period of the ordered phase diverges at the transi- 
tion [23. Thus, the arguments^] that fluctuations 
will generally drive an otherwise continuous freez- 
ing transition first order do not apply; the contin- 
uous character of this transition is robust. 

• 2) There is, of course, some coupling between the 
translational motion of the crystalline order in 
neighboring stripes, so at mean-field level the crys- 
talline order will be locked from stripe to stripe. 
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Consequently, the stripe ground-state breaks trans- 
lation symmetry not only in the direction perpen- 
dicular to the stripes, but along the stripe direction 
as well. However, near the transition, where i<1, 
the spacing between stripes is large compared to 
Lq, so this coupling is exponentially small; conse- 
quently, this locking can be neglected for all prac- 
tical purposes. Therefore, this phase should op- 
erationally be classified as an electron smectic 17], 
in which translation symmetry is unbroken along 
the stripe direction. (There remains the interest- 
ing academic question of principle whether or not 
quantum fluctuations are able to truly stabilize this 
smectic phase at T = - this is closely related to 
the issue of whether "floatin g ph ases" are stable in 
quasi-lD electronic svstems |l7ll2^ '). 

• 3) Near x — 1/2, the stripes of Wigner crystal 
and the intervening Fermi liquid are comparable in 
width. As x — > 1, the system is better thought of as 
stripes of Fermi liquid separated by broad regions 
of Wigner crystal. At some point, the crystalline 
order becomes so rigid that the coupling across the 
liquid stripes is no longer negligible. At this point, 
the striped state is fully crystalline, in the sense 
that translation symmetry is broken in both di- 
rections, and the structure factor contains Bragg 
peaks. However, this phase is still qualitatively 
distinct from the Wigner crystal. Since generally 
speaking the Fermi wave vector is unrelated to the 
Bragg vectors of the Wigner crystal, the liquid in 
the stripe "rivers" can still conduct current in the 
stripe direction. For want of a better name, we 
christen this state a striped "conducting crystall." 
(See Fig.2). 

• 4) At x = 1, the transition from the conducting to 
Wigner crystal mirrors the smectic to Fermi liquid 
transition; it is also a Lifshitz transition at which 
the period of the stripe order diverges. 



B. Bubble Phases 

So far, this analysis ignores the possibility of bubble 
phases. Whether or not there is a regime in which the 
lowest energy mean-field state is a bubble phase depends, 
as we mentioned before, on the degree of anisotropy of 
the microscopic surface tension, ao(9). It may happen, 
due to the anisotropy of the Wigner crystal, that ao(6) is 
sufficiently anisotropic that bubble phases never intrude 
upon the phase diagram. It is also possible to force the 
issue by artificially enhancing the anisotropy of <Jq{0). 
This can be done by explicitly breaking the rotational 
symmetry of the 2DEG, for instance by applying an in- 
plane magnetic field or by using a sufficiently anisotropic 
surface in the construction of the MOSFET. In this case, 
no more need be said. 



However, the Wigner crystal is generally thought to be 
triangular. In this case the surface energy is sufficiently 
isotropic that for x near or 1, there will be a range 
of x in which bubble phases have lower energy than the 
stripe phase; for x near 0, such a phase consists of far 
separated crystallites in a metallic sea, while for x near 
1, it is far separated bubbles of fluid in a Wigner crys- 
talline host. We will call these phases Bubble Crystals I 
and II, respectively. (See Fig. 2d) As x — > or x — > 1, the 
period of the bubble crystals diverge, leading at mean- 
field level to another Lifshitz transition, much as in the 
stripe case. (However, fluctuation effects are much differ- 
ent near these transition in the bubble and stripe cases, 
as we will discuss in the next sections.) 

However, this is not the end of the story. The stripe 
phase is always the lower energy one near x — 1/2, so 
if a bubble phase occurs for small x, there must be a 
critical value of x — x c at which the energy of the bub- 
ble and stripe phases cross, seemingly implying a first 
order transition. Since we have proven in general that 
first order transitions are forbidden, this first order tran- 
sition, too, must be replaced by a regime of intermediate 
phases consisting of a mixture of bubble and stripe phases 
|14j. Now, however, because of the large anisotropy of 
the stripe phase, the surface tension between these two 
phases must be highly anisotropic. Thus, this intermedi- 
ate phase will most probably be of the form of alternating 
mega-stripes of bubble and stripe phase regions. These 
regions are shown in Fig. 2d by hatched boxes. 



II. THERMAL EFFECTS. 
A. The Pomaranchuk Effect 

The most dramatic effect of finite T is its effect on the 
balance between the liquid and Wigner-crystalline phases 
- the fraction of the Wigner crystal phase grows as the 
temperature increases |14)! This phenomenon is similar 
to the Pomeranchuk effect in He 3 and has the same ori- 
gin: the spin entropy of the crystal phase is substantially 
larger than that of the liquid state. The same considera- 
tions lead, as well, to the conclusion that the crystalline 
phase is preferred relative to the liquid in the presence of 
an in-plane magnetic field, H\\ . 

Due to the Pauli exclusion principle, an effective ex- 
change energy of order the Fermi energy, E F ~ h 2 n/2m* 
quenches the spin entropy in the liquid phase. In con- 
trast, the exchange!^ energy in the Wigner crystal is 
exponentially small, J cx exp[— a^/r^] where a is a num- 
ber of order 1. For example an estimate made in [2(j 
yields J ~ 10~ 7 Ry* where Ry* = e 4 m*/2h 2 e is the ef- 
fective Rydberg. Thus a combination of the quantum 
character of the liquid and the smallness of exchange 
processes in the solid imply that the solid phase is sta- 
bilized by non-zero T relative to the liquid phase - for n 
near n c , the 2DEG freezes upon heating! In the present 
context, this means that for fixed n, the relative frac- 
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FIG. 2: Schematic representation of the sequence of inter- 
mediate states as the 2DEG evolves from the Fermi liquid 
(FL) to the Wigner crystal (WC). In a, b, and c, we assume 
that cr{0) is sufficiently anisotropic that bubble phases are su- 
pressed. a) The mean-field phase diagram. Under appropriate 
circumstances, this also represents the true sequence of phase 
transitions at T — 0. b) The phase diagram at non-zero tem- 
perature with a rotationally invariant Hamiltonian. c) The 
phase diagram at non-zero temperature with a preferred axis, 
for instance due to an in-field magnetic field. Vertical lines 
represent phase transitions and wavy lines crossovers. Phases 
with power-law order are named in italics. The smectic phase 
in a) is weakly unstable to crystallization at mean-field level, 
but may be stabilized by quantum fluctuations, d) The phase 
diagram including bubble phases at T — in the presence of 
quantum fluctuations. The hatched areas correspond to the 
sequence of transitions involving mega-stripes of bubble and 
stripe phases discussed in the text. 



tion of Wigner crystalline regions increases with increas- 
ing T or Hu . A simple estimate of the magnitude of 
this effect can be made for the range of temperatures 
J < T < E* F = Ry* /ir{r s ) 2 and h/j, B H\\ < E* F , where 
the entropy of the liquid is negligible, as are the subtleties 
of the ground-state magnetic structure of the Wigner 
crystal. In this case, 

f(n-,T,H\\) «/(n_,0,0) (8) 
-k B Tn^\n [2cosh(hn B H {l /k B T)] 

and f(n+,T,H\\) w /(n+,0,0). The fact that tempera- 
ture and magnetic field stabilize the Wigner crystal in 
qualitatively similar fashion is one of the striking as- 
pects of this relation: For T ^> hfi B H\\, f(n-,T,H\\) — 
f(n-, 0,0) -/c B Tn_ln[2] while for T < h/i B H h 
/(n_,T, B) - /(n_, 0, 0) w -ft/i B #||n_. 

Of course, at high enough temperatures, all tenden- 
cies to ordered states are suppressed. This occurs above 
the characteristic temperature at which the Wigner crys- 
tal melts. In the limit of very large r s , this occurs at 
the classical melting temperature of the Wigner crys- 
tal, which has been estimated in accurate numerical 



experiments plj to be 

T melt = A(e 2 /e)V^ = 2ARy*/r s (9) 

where A = 1/125 [1± 0.04]. However, at smaller r s , where 
Ef of the competing fluid phase is larger than the puta- 
tive classical melting temperature, the implied reduction 
of the entropy of the fluid state means that the melt- 
ing temperature is set, by T me / t oc Ep. Far from the 
Lifshitz points, the melting temperatures of the various 
microemulsion phases are determined by these same con- 
siderations, and are of similar magnitude. Here, the frac- 
tion of the system that is crystalline is a non-monotonic 
function of T, first increasing and then dropping to zero 
at T me i t . Near the Lifshitz points, more delicate consid- 
erations determine the melting point. 

B. Thermal fluctuations in the stripe phases 

Let us now consider the role of thermal fluctuations 
on the stripe phases. We distinguish two cases: 1) If the 
Hamiltonian is rotationally invariant, then the smectic 
phase is unstable at any non-zero temperature to the 
proliferation of dislocations. Thus, 

the mean-field smectic phase is replaced by a nematic 
phase, which, in keeping with the Mermin- Wagner the- 
orem, does not actually break rotational symmetry, but 
rather has power-law orientational order. A free disloca- 
tion has a logarithmically divergent energy in both the 
Wigner and conducting crystal phases, so they are ro- 
bust against thermal fluctuations at low temperatures, 
although with power-law rather than long-range crys- 
talline order. The resulting phase diagram is shown 
schematically in Fig. 2b. 2) If, however, the Hamilto- 
nian has a preferred axis, for instance if we consider the 
2DEG in the presence of an in-plane magnetic field, the 
effects of low temperature thermal fluctuations are much 
less severe. Here, the smectic and both crystalline phases 
remain well defined at non-zero T, although again with 
power law spatial correlations rather than with true long- 
range order, as shown schematically in Fig. 2c. 

Because first order transitions are forbidden, the tran- 
sition between the isotropic fluid and the nematic phase 
must be of the Beresinskii-Kosterlitz-Thouless (BKT) 
type. Near the mean field Lifshitz point we can estimate 
this transition temperature as follows: The distance be- 
tween stripes is large so the stripes of minority phase 
evaporate when the energy to break off a piece, ~ ctqLo, 
is less than the the configurational entropy of a state 
where rare droplets of the minority phase are distributed 
randomly. Equating these two free energies leads to the 
estimate 

T c |m[x(l-x)]| ~ ao L . (10) 

In the presence of an in-plane magnetic field, there is no 
sharply defined nematic phase, since rotational symme- 
try is explicitly broken. However, by the same token, free 
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dislocations in the smectic state have a logarithmically 
divergent energy, so a power-law smectic phase exists at 
non-zero T. With increasing temperature, the smectic to 
liquid phase transition is also of the BKT type. Indeed, 
so long as the symmetry breaking term in the Hamilto- 
nian is small, the transition temperature is roughly the 
same as in Eq. El above. 



C. Thermal fluctuations in the bubble phases 

We now consider the effect of thermal fluctuations on 
bubble phases. As can be seen from Eq. 1, the interaction 
energy between far separated bubbles decreases at large 
r as 

Vbubble ~ eriLo/rj], (11) 

where Lq and ro are the radius of and distance be- 
tween bubbles. Thus, where the bubbles are far sepa- 
rated, because of the screening by the ground-plane the 
BKT melting temperature will tend to rapidly to zero, 
Tbkt oc [(1 — x)x] 3 / 2 , as the spacing between bubbles 
increases. The result is that, near the mean filed Lifshitz 
point the bubble phase is always melted by the thermal 
fluctuations. On the other hand, at smaller r the bub- 
ble phase survives thermal fluctuations in the usual sense 
that the correlations of bubble positions exhibit power- 
law decay. 

The nature of the transition between the bubble phase 
and the uniform phase is not, presently, settled. Of 
course, a direct first order transition is forbidden. One 
possibility is that there is a sequence of two BKT tran- 
sitions, as in the Halperin-Nelson theory[2^ of melt- 
ing, with an intermediate hexatic phase. Alternatively, 
there may be a further set of hierarchical microemulsion 
phases. 



The quantum nature of the system near the Lifshitz 
points is determined by the quantum nature of the in- 
terface between the crystal and the liquid - a problem 
which itself is still unsolved. This interface may be quan- 
tum rough or quantum smooth. If it is smooth, the Lif- 
shitz transition from the uniform fluid to smectic phase 
is not fundamentally affected by quantum fluctuations, 
provided the width of the stripes is large enough. How- 
ever, if an isolated interface is rough, the stripe order 
in the vicinity of the mean- field Lifshitz point is quan- 
tum melted; in this case, for the rotationally invariant 
system, the proscription against first order transitions 
implies that there must be an intermediate zero temper- 
ature nematic phase between the isotropic and the stripe 
ordered phases. It was recently shown 18] that a nematic 
Fermi fluid is necessarily a non-Fermi liquid in the sense 
that quasiparticles are not well defined elementary ex- 
citations. We believe that, depending on microscopic 
details and on the the value of nd 2 , both scenarios are 
possible. 

It is worth mentioning why the quantum nature of the 
crystal-liquid interface is so subtle. Consider the mo- 
tion of a step in the interface. Quantum-mechanically, 
an isolated s tep might be expected to propagate along 
the interface [23. Because the steps interact by a short- 
range dipolar interaction, the steps should then form a 
delocalized ID quantum liquid along the interface. How- 
ever, because of the density mismatch between the solid 
and liquid, the situation is more complicated. Motion of 
the step requires a flux of mass into the liquid of a mag- 
nitude proportional to the density difference between the 
solid and liquid and to the step's velocity. In a Fermi liq- 
uid this flux of mass is carried by quasiparticles, making 
the step motion highly dissipative. Thus, characterizing 
the interface involves interesting, but as far as we know 
unsolved issues in dissipative quantum mechanics. 

B. Bubble Phases 



III. QUANTUM FLUCTUATIONS 
A. Stripe Phases 

So long as nL\ 1 (nd 2 ^ 1), the stripes are many 
electrons wide, so quantum fluctuations of their positions 
are intrinsically small; 1 / nL\ is a small parameter in the 
problem, which permits an asymptotically exact treat- 
ment of quantum fluctuation effects. At zero temper- 
ature, the conducting-crystal phase is clearly stable in 
the presence of small quantum fluctuations, although, as 
mentioned previously, the jury is still out on whether the 
smectic phase is unstable to crvstallization |l7ll22l ]. Only 
where the stripe width is of order of the interelectron 
distance ( i.e., when nL\ ~ 1), quantum fluctuations be- 
come very significant. This applies to the hatched region 
in Fig. ^ where the quantum properties of the system 
are still uncertain. 



In contrast to stripe phases, quantum fluctuations al- 
ways melt the bubble phases when the bubbles are suf- 
ficiently dilute. To see this, we can estimate the charac- 
teristic potential energy of a bubble crystal as in Eq. 1111 
and can make a corresponding dimensional estimate of 
the bubble kinetic energy Kbubbie ~ h 2 /r^m* where m* 
is the bubble effective mass. Therefore the ratio of these 
energies is 

Vbubble ^ f (TiL%m* \ 1_ 
Kbubbie V h 2 J r 

vanishes as r$ — > 00. This analysis fleshes out the same 
argument mentioned in the introduction that leads to the 
conclusion that there is no stable Wigner crystal phase at 
small d. However, whereas in that case, the proportional- 
ity constant is a* B , in the present case the same constant 
is paramctrically large, both due to the explicit factors of 
Lq and due to the fact that m* increases with increasing 
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Lq (in a way that depends on whether the interface is 
quantum rough or smooth). The result is that, for large 
nd 2 , the regime in which the bubble crystal is quantum 
incited is extremely small. However, if rid 2 ~ 1 quantum 
melting is a significant phenomenon. 

The character of the bubble liquid phase is different 
depending on the character of the minority phase. When 
the majority phase is Wigner-crystalline with dilute in- 
clusions of liquid, the melting of the bubble crystal re- 
sults in a type of "conducting crystal" 0] . In this state, 
crystalline long-range order coexists with fluid-like con- 
ductivity, but in this case the conductivity is associated 
with the motion of the bubbles themselves. Phenomeno- 
logically this state is similar to the "supersolid" phase 
which has been discussed[2|| in the context of He 4 . In 
both cases the number of electrons (or aroms) is not equal 
to the number of the crystalline sites. The differ- 

ence is that unlike the case of He 4 where vacancies are 
bosons, in our case the statistics of the droplets is not 
known, and hence the liquid state may not be a super- 
fluid. Therefore, we refer to this state as a "Conducting 
Crystal II" in Fig. 2d to distinguish it from the highly 
anisotropic conducting crystal (See Fig. 2a, d) which orig- 
inates from the existence of stripes. 

When the Fermi liquid is the majority phase, with a 
fluid of "icebergs" floating in it, no spatial symmetries 
need be broken. However, elementary excitation spec- 
trum is likely to be different from that of a conventional 
Fermi liquid. 

Since the majority phase already brakes rotational 
symmetry, the bubble crystallization transition which 
transforms the system from the conducting crystal to the 
insulating bubble crystal phase can be a simple continu- 
ous transition. However, the freezing of the icebergs into 
a triangular crystal of Wigner-crystalline bubbles is more 
problematic. As with the thermal transition, there may 
be a two-step freezing transition, with an intermediate 
quantum hexatic phase. IS], or another hierarchy of mi- 
croemulsion phases. The sequence of the phases at T = 
is shown in Fig. 2d. 



IV. EXPERIMENTAL CONSEQUENCES 

Obviously, there are many experimental consequences 
of the existence of intermediate phases, of which we here 
list only a few. It should be kept in mind that macro- 
scopic spatial symmetry breaking, the sort which pre- 
cisely characterizes the various phases we have discussed, 
does not truly occur in 2D in the presence of quenched 
disorder. This complicates the actual observation of var- 
ious phenomena. 

The majority of industrially produced Si MOSFET's 
have gates relatively close to the 2DEG, d <c d c , so the 
electron liquid is weakly interacting at all n. However, 
a small number of high mobility Si MOSFET's (For a 
review, see |30|.) and p-type of GaAs double layers 
with large d ~ lOOOA have been studied in the past few 



years, and found to exhibit transport anomalies that have 
been interpreted as evidence for an unexpected metal- 
insulator transition. While these devices certainly are 
not ideal, in the sense that they have non-zero quenched 
disorder, we would like to propose that a natural explana- 
tion of these phenomena is that they reflect the existence 
in the zero disorder limit of the electronic microemulsion 
phases identified in the present paper. 

One robust consequence of two-phase coexistence is 
that the conductivity is a decreasing function of the vol- 
ume fraction of Wigner crystal. This volume fraction, 
in turn, is strongly temperature and magnetic field de- 
pendent due to the Pomeranchuk effect, as explained 
above. As a result, the fraction of crystal grows with 
temperature and magnetic field, leading to a correspond- 
ing increase of the resistivity. As has been pointed out 
previously this basic physics may underly the trans- 
port anomalies observed in large d Si MOSFETs. In par- 
ticular, it offers a candidate explanation of the anoma- 
lous metallic (dp/dT > 0) temperature dependence and 
large positive magneto-resistance observed in these sys- 
tems despite the fact that the resistivity, itself, exceeds 
the Ioffe-Regel limit. (Ideally, one might want to explore 
the scaling relation between the temperature and mag- 
netic field dependence of the resistivity implied by the 
thermodynamic relation in Eq. [5J) 

Of course, each new phase has different patterns of 
spatial symmetry breaking, and hence has new collective 
modes and modified hydrodynamics. Even when the ef- 
fects of quenched disorder or thermal fluctuations restore 
the symmetry at macroscopic distances, the existence of 
these collective modes can have readily detectable conse- 
quences for the dynamical responses of the system. Small 
explicit symmetry breaking fields can be used to over- 
come the destructive effects of quenched disorder and 
reveal the true tendency to symmetry breaking. For in- 
stance, an in-plane magnetic field explicitly breaks rota- 
tional symmetry; where some form of stripe or nematic 
phase exists in the absence of quenched disorder, the 
small symmetry breaking produced by such a field can 
give rise to a large resistivity anisotropy, as has been seen 
for the analogous states in quantum Hall devices |24j. 



V. EXTENSIONS 

We end with some speculative observations concerning 
intermediate phases of the 2DEG. 

Spin physics: Other than the Pomeranchuk effect, we 
have largely ignored the physics of the electron spins. 
The exchange interactions in the Wigner crystal phase 
are generally found to be very small|25|], and so are only 
important at very low temperatures. At T = 0, how- 
ever, the fact that the magnetic Hamiltonian is highly 
frustrated and may have important multispin ring ex- 
change interactions, can lead to a variety of possible mag- 
netic phases, and this complexity could be inherited, to 
some degree, by the intermediate phases discussed here. 
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Moreover, at a liquid-crystalline interface, the quantum 
dynamics of the interface itself (mentioned above) can 
produce effective exchange interactions, likely with much 
larger energy scales than in the bulk Wigner crystal. 
There is thus the very real possibility that the magnetic 
structure of the interfaces is very rich, and characterized 
by substantial energy scales. 

Superconductivity: The parallels between the 2DEG 
in a MOSFET and Coulomb frustrated phase separation 
in a doped Mott insulator naturally lead one to spec- 
ulate concerning the possibility of superconductivity in 
the present system. In the bubble related conducting 
crystal phase, each bubble has a fixed number of elec- 
trons; when that number is even, the bubbles are likely 
bosonic and a supersolid phase with low superfluid den- 
sity is possible [lj,|28|]. In the hatched region of the phase 
diagram, where quantum effects are most severe, a more 
robust mechanism is possible, based on the "spin-gap 
proximity effect 26]": Small clusters of Wigner crystal 
(be they stripe or bubble like) will often have a spin-gap. 
(Near the cluster edge, this gap may be larger than in 
the bulk.) Where this gap is large enough, it suppresses 
single-particle exchange between the crystal clusters and 
the surrounding Fermi fluid, but pair-exchange is still 
permitted. When this dominates, it induces global super- 
conductivity by a process analogous to the conventional 



proximity effect. 

Double Layers: In a double layer system, with two 
nearby 2DEG's, the two layers screen each other in much 
the same way as the metal layer screens the 2DEG in 
a MOSFET. However, here the types of phases, and the 
available experiments are still richer. One particularly in- 
teresting point is that the conductivity measured in drag 
can explore the nature of the interlayer screening. The 
presence of a crystalline component of an electron fluid 
has the potential to greatly increase the drag conductiv- 
ity relative to a Fermi liquid; in particular, whereas the 
drag conductivity vanishes as T — > in a Fermi liquid, 
we believe it can approach a non-zero constant in some 
of the intermediate phases we have explored. 

Other Applications: The present ideas are pretty 
clearly applicable in a host of additional physical con- 
texts. What is needed is short-range tendency to phase 
separation, i.e. a concave local free energy, opposed by 
dipolar forces. Under appropriate circumstances, this 
situation may pertain in the 2DEG at higher densities, 
r s < r c , and it certainly applies in various regimes to the 
physics of lipid films and planar ferromagnets. 
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